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ABSTRACT 

We review the evidence for substructures from the anomalous flux ratios in 
gravitational lenses. Using high-resolution numerical simulations, we show that 
at typical image positions, the fraction of surface mass density in substructures is 
< 0.5% with mass above 10~^ virial masses in the "concordance" ACDM cosmol- 
ogy. Substructures outside the virial radius (but projected at typical lens image 
positions) only increase the fraction moderately. Several effects, in particular 
baryonic settling and the requirement of compactness, may further decrease the 
predictions by a factor of few. The predicted fraction with appropriate prop- 
erties thus appears to be lower than that required by lensing, although both 
are still uncertain. More speculative substructures such as massive black holes 
(M ~ 10^ — IO^Mq) in the halo may offer viable alternatives. 

Subject headings: gravitational lensing - cosmology: theory - dark matter - 
galaxies: structure, evolution 



1. INTRODUCTION 

Gravitational lenses on arcsecond scales provide a unique sample to probe the mass 
distribution in the lensing galaxies at intermediate redshift {z ~ 0.5 — 1). Image positions 
in most lenses can be fitted adequately using simple smooth galaxy mass models. But 
observed flux ratios are more difficult to match (e.g., Kochanek 1991). The discrepancy 
between the predicted and observed flux ratios is commonly referred to as the "anomalous 
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flux ratio problem." The most apparent cases are found in quadrupole lenses where we 
observe a close pair or a close triple of images. Here we know that the lensed source is close 
to either a fold or a cusp caustic. The asymptotic magnification behavior in such cases is 
well understood - a close pair must have equal brightness, while for a close triple, the flux 
of the middle image should be equal to the total fluxes of the two outer images. Virtually 
all the observed pairs and triples disagree with these relations (§2). This has been argued 
as evidence for substructures on the scale of the separations of the images (a few tenths of 
arcseconds, e.g., Mao & Schneider 1998; Metcalf & Zhao 2002). Another piece of evidence 
for substructures is the fact that saddle images are preferentially dimmed compared to model 
predictions (Kochanek Sz Dalai 2003). This is expected from miUi-lensing by substructures 
(Keeton 2003) or microlensing by stars (Schechter & Wambsganss 2002)^; such a preferential 
de-magnification of saddle images is not expected from other propagational effects. 

The Cold Dark Matter (CDM) structure formation model predicts the existence of just 
such substructures from both semi-analytical studies and numerical simulations (e.g., Kauff- 
mannetal 1993; Klypin ct al. 1999; Moore et al. 1999; Ghigna ct al. 2000). About 5-10% of 
the mass is predicted to be in substructures with a mass spectrum of n{M)dM ~ M'^-^dM. 
Intriguingly, the predicted number of subhaloes in CDM exceeds the observed number of 
luminous satellite galaxies in a Milky- Way type galaxy (e.g., Klypin et al. 1999; Moore et 
al. 1999; see also Stoehr et al. 2002). One solution for this disagreement may be that some 
substructures, especially those of lowest mass, are dark. If this is true, then gravitational 
lensing may be the best way to detect them. In this paper, we examine the required amount 
of substructures in gravitational lenses (§2) and compare it with the predictions from CDM 
(§3). Finally, in §4, we discuss several effects that affect the predictions of numerical simula- 
tions. Throughout this paper, we adopt the "concordance" ACDM cosmology (e.g., Ostriker 
& Steinhardt 1995; Spcrgel et al. 2003 and references therein), with a density parameter 
^m,o = 0.3, a cosmologically constant ^a,o = 0.7, a baryon density parameter f2b = 0.024/?."^, 
and we take the power-spectrum normalization ag = 0.9. We write the Hubble constant as 
h = i/o/(100kms-iMpc-i) with h = 0.7. 

2. Required Amount of Substructures From Gravitational Lenses 

A number of papers have discussed evidence for substructures from gravitational lenses; 
most concentrated on the anomalous flux ratios (e.g., Mao & Schneider 1998; Metcalf & 



^We refer to lensing by substructures and lensing by stars as milli-lensing and microlensing respectively 
because the angular scales involved are ~ mas and fias in these two cases. 
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Madau 2001; Chiba 2002; Bradac et al. 2002; Dalai & Kochanek 2003; Keeton 2003), while 
several papers discussed astrometric signatures, such as bent jets in B1152+199 (Metcalf 
2002), and unusual VLBI structures for MG2016+112 (Kochanek & Dalai 2002; see also 
Koopmans et al. 2002) and B0128+437 (Phillips et al. 2000; Biggs et al. 2003). It is 
illustrative to see the issues using the largest homogeneous lens survey - the Cosmic Lens 
All-Sky Survey (CLASS, Myers et al. 2003 and Browne et al. 2003). This radio survey has 
well defined selection criteria and does not suffer from the effect of dust extinction. Radio 
lenses are also not substantially affected by stellar microlensing (see Koopmans et al. 2003 
and references therein). 

In total, there are 22 new lenses discovered in CLASS. Among these, 7 are simple 
quadrupole lenses, including 5 close pairs (B0128+437, Phillips et al. 2000; MG0414+0534, 
Hewitt et al. 1992; B0712+472, Jackson et al. 1998; B1608+656, Myers et al. 1995; 
B1555+375, Marlow et al. 1999), and two close triples (B2045+265, Fassnacht et al. 1999; 
B1422+231, Patnaik et al. 1992). For the five pairs listed above, the observed fiux ratios are 
0.56 (5GHz), 0.88 (15GHz), 0.75 (15GHz), 0.51 (8.4GHz), 0.56(15GHz); for each lens, the 
highest frequency where the pair has been observed is shown in the bracket. For the two close 
triples, the ratios of the flux of the middle image to the total flux of the two outer images 
are, 0.32 (14.9GHz) and 0.70 (15GHz, Patnaik & Narasimha 2001) respectively. The typical 
errors on the fiux ratios are ~ 1 — 2%. Except the case of MG0414-I-0534, the observed 
values are all different from the asymptotic value (unity). 

However, not all these flux ratios are "anomalous" because some of these systems may 
not yet reach the asymptotic regime. Also some radio lenses may have been affected by 
scattering by free electrons along the line of sight. In addition, most lens modelers adopt 
isothermal ellipsoidal models or variants as an approximation for the lensing potentials. 
Models with more complex radial and angular structures can usually better reproduce the 
observed fiux ratios (e.g., Evans & Witt 2003; however, see Kochanek & Dalai 2003). Obser- 
vationally, it is clear that some lens systems are complex. For example, B1608+656 has two 
lensing galaxies and a model which accounts for this fact can reproduce the observed flux 
ratio (Koopmans & Fassnacht 1999). For B0712+472, a singular isothermal ellipsoid (SIE) 
model cannot reproduce the flux ratio, but a foreground group of galaxies was subsequently 
found along the fine of sight of the primary lensing galaxy (Fassnacht & Lubin 2002), so 
a more realistic model may match the fiux ratio. The close pair fiux ratio in B1555-I-575 
cannot be reproduced by an SIE model (Meyers et al. 1995), but a more complex model 
(e.g., with an additional shear) may be able to explain it. For B0128-I-437, an SIE model 
cannot fit the fiux ratio, but the observed fiux ratio can be reproduced with an additional 
shear. However, in this case, the relative orientations of jets in this system appear difficult 
to match with a smooth model (Biggs et al. 2003), so substructures may yet be called for. 
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For the two close triple lenses, the flux ratios appear difiicult to match with smooth models 
(Mao & Schneider 1998; Fassnacht et al. 1999). 

Mao & Schneider (1998) showed that substructures preferentially affect the flux ratios 
of highly magnifled images. In order to reproduce the flux ratios in B1422+231, they found 
that the required perturbation in the dimensionless surface mass density, 5k/ is of a few 

percent and roughly corresponds to a physical surface density of a few tens MqPc"^. Dalai 
& Kochanek (2002) performed a statistical study of seven radio lenses using Monte Carlo 
simulations, six of which show anomalous flux ratios. They find that the fraction of mass 
required in substructures should be in the range of /sub = 0.6% — 7% (90% confidence limit) 
with a best fit of /sub ~ 2%. More recently, Metcalf et al. (2003) applied the method of 
Moustakas & Metcalf (2003) to the quadrupole system 2237+0305. They conclude that, in 
order to match the observed flux ratios in the radio, infrared and narrow and broad emission 
fines, 4%-7% of the surface mass density (95% confidence limit) must be in substructures 
with mass between lO^Mp, - 10^M„. 



3. Predicted Substructures in Numerical Simulations 

We use the high resolution halo simulations of Jing & Suto (2002; 2000) to constrain the 
fraction of mass in substructures; similar results (but with larger error bars) are found using 
the simulation data obtained with the tree-particle- mesh code of Bode & Ostriker (2003). 
Twelve halos were selected from a cosmological simulation of box of 100/i~^ Mpc, with four 
each at galactic, group, and cluster masses, respectively. These simulations are evolved with 
a nested-grid PPPM code which was designed to simulate high-resolution halos. The force 
resolution is typically 0.004rvii., where ryij. is the virial radius. At the end of each simulation, 
about (0.5 — 1) X 10^ particles are within of each halo (see Jing & Suto 2000 for details). 

We adopt the SUBFIND routine of Springel et al. (2001) to find disjoint self-bound 
subhalos within a parent halo. All subhalos with more than 10 particles are included in our 
analysis. Our numerical resolution implies that we can only identify subhaloes with mass 
larger than about 10^^ of the parent halo mass; the most massive subhalo has about 10% of 
the parent halo mass. For each halo, we make 30 random projections and calculate the total 
mass within different annuli; the annuli are equally spaced in log R/r^u from —2.2 to with a 
step size of 0.2, where it! is the projected radius and the spherical radius is r. The total mass 
in substructures within an annulus is calculated by summing up all the mass of subhalos 
whose centers fall in it. Dividing the total substructure mass by the total mass within the 
annulus yields the fraction of the projected surface density in substructures, /sub- The mean 
value and variance of /sub are found from the 12 haloes and the 30 random projections. In 
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the upper left panel of Fig. 1, we show /sub as a function of R/r^ir- The mean fraction can 
be approximated as f^uh ~ 0.25(i?/rvir). The scatter around the mean among haloes is quite 
large. At i? ~ ryir, the scatter is about 40%; at smaller R the scatter is larger due to fewer 
subhaloes along the line of sight. As lensing concerns only the projected surface density, 
wc also checked whether substructures outside the spherical virial radius can contribute to 
the surface density at a given R (Fig. 1, top right). The increase is about a factor of two 
R ^ Tvir (cf. Klypin et al. 2001), but more modest (< 1.3) at typical image positions 
{R ~ O.Olrvir — 0.03rvir), consistent with the analytical estimate by Chen et al. (2003). 

Almost all subhaloes in the simulation have a spherical radius r > O.lrvir- Most sub- 
structures at typical image positions {R ~ O.Olrvir — 0.03rvir) only appear along the line of 
sight due to projection. These subhaloes, especially those at r ~ Tyir, are extended, thus 
not all the bound mass within the subhaloes can efficiently cause the flux anomaly (see §4) . 
To illustrate this effect, /sub is recalculated by including only the mass within a spherical 
radius of 0.025rvir for each subhalo. Due to this compactness requirement, at i? = 0.03rvir, 
the average /sub value (for substructures with r < Tyir) is lowered from 0.8% to 0.6% (see the 
lower panels in Fig. 1). We return to this compactness issue in the discussion. 

To compare the predicted fraction with lensing requirement we need an estimate of the 
virial radius for the lensing galaxies. If the lens and source redshifts arc known, then the 
separation between images (oc o"^) can be used to determine the velocity dispersion (a) and 
the approximate halo circular velocity {V^ ~ v^cr), which in turn allows us to determine the 
halo virial radius (e.g., cqs. 2-3 in Mo, Mao & White 1998). The lens and source redshifts 
are known for MG0414+0534, B0712+472, B1608+656, B2045+265 and B1555+375. For 
these systems, the image positions in units of the virial radius are, 0.027, 0.012, 0.018, 0.017 
and 0.012, respectively. From Fig. 1 we can then infer the fraction of mass in substructures 
is only around 0.2-0.8%, with a large scatter among different haloes. 

So far we have only considered the fraction of dark matter surface density in substruc- 
tures. However, the inner parts of galaxies are likely dominated by baryons. To address this 
issue, we consider the coohng and settling of baryons in an NFW halo (Navarro, Frenk & 
White 1997), similar to the procedure used by Keeton (2001). Initially the baryons and dark 
matter follow the same profile. The baryons then cool and condense into the center to form 
a stellar component with a de Vaucouleurs profile while the dark matter responds to the 
baryonic settling adiabatically. The model is described by the concentration parameter c in 
the NFW profile, the effective radius (^cfr) in units of the halo virial radius, and the ratio of 
the stellar mass to the total mass, /cooi (or equivalently, a mass-to-light ratio). 

For the five sources with known lens and source redshifts (MG0414+0534, B0712-I-472, 
B1608-I-656, B2045-I-265 and B1555-I-375), the effective radii have been determined using 
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HST photometry (Kochanek et al. 2000). In units of the virial radius they are, respectively, 
0.018, 0.006, 0.01, 0.006 and 0.005. We take c = 10 (appropriate for a galactic-sized halo) 
and then carry out the procedure described above. For the five systems, we find that the 
stellar component contributes about 20%-50% of the projected surface density at the image 
positions (i? ^ l-Srcff — 3rcfr) for /cool ranging from 0.05 to 0.16 (~ f2b/^o)- The fraction of 
surface mass density in dark matter is consistent with that required by stellar microlensing 
(~ 70 — 90% at l.Sreflf, Schechter & Wambsganss 2003), but at odds with the recent claim 
of Romanowsky et al. (2003). The effect of baryons hence reduces the the mass fraction in 
substructures in typical lensing systems to < 0.5%. 



4. Discussion 

We have reviewed the evidence for substructures from close pairs and triples in quadru- 
ple lenses. As emphasized by Kochanek & Dalai (2003), the fact that saddle images are 
frequently dimmer than expected is difficult to accommodate by other means. Quantita- 
tively, the anomalous flux ratios in lenses appear to require a few percent of the surface mass 
density in substructures at typical image positions (Dalai & Kochanek 2002; Metcalf et al. 
2003). The required fraction is higher than that provided by globular clusters and luminous 
satellite galaxies (Mao & Schneider 1998; Chiba 2002) and it also appears to be higher than 
the predicted values (/sub ^ 0.5%) from the ACDM cosmology at typical image positions. 
However, at present it is unclear how serious the discrepancy is because of uncertainties in 
both observations and theoretical predictions. 

There are a number of issues that need to be understood better in current numerical 
simulations. Even the basic question of the identification of substructures needs to be ex- 
plored further. The SUBFIND algorithm we adopted only identifies bound subhaloes, however, 
tidal streams from disrupted systems (for examples in the Milky Way, see, e.g., Ibata et al. 
2001; Yanny et al. 2003) may also contribute to the budget of substructures. Another impor- 
tant issue is whether our results have achieved convergence as a function of the spatial and 
mass resolutions. New simulations are underway to address this issue. Presumably when 
the numerical resolution becomes higher, the inner parts of subhaloes are resolved better 
into higher-density regions that can survive tidal disruptions longer. However, the survival 
of substructures may be linked to another small-scale problem of CDM: the central density 
profiles of low-surface brightness galaxies (usually with circular velocities of < 100 km s"^) 
seem to be too concentrated compared with observed galaxies (e.g., Bolatto et al. 2002; 
Weldrake et al. 2003; see, however, Swaters et al. 2003). Therefore, if we put in observed 
mass profiles, substructures may actually be more easily destroyed by tidal forces due to their 
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lower central concentrations. There is another effect that makes the survival of substructures 
in the central part more difficult. In collisionless numerical simulations, the density profile 
can be approximated by an NFW profile; the density scales as oc out to ~ O.lrvir (~ 25 
kpc) for typical galactic-sized haloes. However, the observed velocity dispersion is nearly 
constant in the inner part implying p oc r^^, i.e., the density in real galaxies rises more 
quickly as the radius decreases. Hence substructures will be disrupted more easily if they 
come close to the center and dynamical frictions dragging them into the center would also 
be larger in realistically simulated galaxies. 

Our simulations resolve substructure masses from IO^Mq to IO^^Mq for a lO^^Af^ parent 
halo. However, according to Metcalf et al. (2003), the required substructure mass is in the 
range of 10^ Mq — IO^Mq in the case of 2237+0305. If the mass spectrum of substructures 
n{M)dM oc M'^-^dM extends all the way down to 10^ Mq, one can estimate that the fraction 
of mass in substructures with lO^M© < M < lO^M© is about a factor of 3 smaller than that 
in substructures with IQ^Mq < M < M^^Mq, i.e., the mass fraction in substructures in the 
range required by Metcalf et al. (2003) will be even lower than our predicted value. This 
conclusion depends on the mass spectrum. However, for any power-law spectrum with a 
slope shallower than —2, the mass fraction in substructures with IO^Mq < M < IO^Mq is 
likely smaller than that in higher-mass substructures identified in numerical simulations. 

There is another effect that may reduce the utilizable mass in substructures even fur- 
ther: the need for compactness, an issue we already touched upon in §2 (see Fig. 1). In 
order to affect the ffux ratios significantly, the physical size of substructures must be suffi- 
ciently compact (cf. Metcalf et al. 2003). The most efficient substructures should be of the 
same order of the image separation of the pairs or triples. For the five CLASS quadrupole 
lenses that have known lens and source redshifts (see §4), the closest pairs have projected 
separations from 0.6/i-^kpc (for B1422+231) to 2.3/i-^kpc (for MG0414+0534). At redshift 
of 0.5, dark matter haloes with M < 4 x IO^Mq will have Tyir ~ l/i~^kpc, and so they will 
be efficient in causing fiux anomalies if they are located between close pairs or triples. For 
larger masses, the effect is reduced. We estimate the reduction by assuming that only the 
mass enclosed within l/i~^kpc will affect the ffux ratios. We adopt a mass spectrum for dark 
matter haloes of n{M)dM oc M'^-^ for lO^M© < M < lO^^M©. Each halo is described 
by an NFW profile with a concentration parameter given by c ~ ll(M/10^^/i~^)°'^^ (e.g., 
Zhao et al. 2003). Note that the substructures in the central parts are tidally truncated 
and so cannot be fitted well by NFW profiles, but as most substructures have r > O.lryir 
in spherical radius, the effect of tidal truncation may be modest. We find that the mass in 
substructures that can efficiently cause fiux anomalies is reduced by an additional factor of 
five compared with the total mass in all substructures - most substructures which are in the 
outer parts are too extended to cause fiux anomalies efficiently. Our rough estimate shows 
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that the compactness requirement is an issue that needs to be addressed more carefully. 
The two additional effects noted can reduce the likely mass fraction in substructures having 
the required masses and sizes to as low as ~ 0.03%, uncomfortably low compared with the 
observational requirement. 

Progress can be made from both the observational and theoretical fronts to reduce the 
uncertainties. Observationally, in the radio, the effect of scattering can be studied with 
observations at high frequency where it is expected to be unimportant. In the infrared, 
it would be interesting to have more observations with integral field spectroscopy similar 
to that reported by Metcalf ct al. (2003). This method offers a way to separate stellar 
microlensing from substructure milli-lcnsing. More astrometric signatures of substructures 
will be important as well (see §2). Theoretically, higher-resolution simulations are needed 
and are already under-way. If future observations and numerical simulations still indicate 
a discrepancy between lensing requirements and CDM predictions, then alternatives must 
be sought. One possibility is that these substructures are massive black holes with M ~ 
10^ - IO^Mq (Lacey & Ostriker 1985; Xu & Ostriker 1994), which satisfy the mass and 
compactness requirements. We require only a few percent of the surface density in the 
substructures, so the density parameter in these black holes is only /sub^o ~ 0.012(/sub/0.04), 
which is about 30%(/sub/0.04) of the baryon density in the universe. These massive black 
holes will have other observable signatures (e.g., Wambsganss & Paczyhski 1992; Tremaine 
& Ostriker 1999) and can be further tested. 

We thank I. Browne, H. J. Mo, T. York and S. White for helpful discussions, X. Kang 
for identifying subhalos and V. Springel for providing us with his SUBFIND code. YPJ is 
supported in part by NKBRSF (G19990754) and by NSFC. 
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Fig. 1. — Predicted fractions of dark matter surface mass density in substructures, /sub, 
as a function of the projected radius, R, in units of the virial radius, r^ir- The lensed 
images typically have R/r^ir = 1% — 3%. The left panels show /sub as a function of R 
for all substructures within the spherical radius r < rvir while the right panels include 
all substructures with r < 2rvir. In the lower panels we only include the mass within 
0.025rvii. for each subhalo. The solid line in each panel indicates an unweighted least-square 
fit (log /sub = log^ + log-R/r^ir, i.e., /sub = AR/r^ij.). The A values are, 0.25, 0.31, 0.20 and 
0.17, clockwise from the top left to the bottom right panels, respectively. 



